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Abstract 

Recently, the behavior of different epidemic models and their relation both to 
different types of geometries and to some biological models has been revisited . Path 
equations representing the behavior of epidemic models and their corresponding 
deviation vectors are examined. A comparison between paths and their deviation 
vectors in Riemannian and Finslerian Geometries is presented. 

1 Introduction 

In this study, we are going to describe a specific susceptive-infective (Sl)-model of 
epidemics using Antonelli's approach of describing two species in ecology geometri- 
cally [1]. This approach has been successfully used in case of wild rabbit disease as 
well as coral-starfish equilibrium [2- 5]. Also, allometric relationship between net pro- 
duction biomass and the amount of secondary compound plants produce to defend 
against herbivores can be geometrized as well as corresponding interaction condi- 
ments are taken to be constant. The aim of this work is to apply Antonelli's idea 
of imposing geometrical paths defined in metric spaces extracted from stochastic 
allometry space in Riemannian geometry to obtain a better non-linear regression 
formula with geometrical origin [2] . The importance of allometry space, due to 
Antonelli's approach, is to obtain a space having a positive definite metric with 
negative curvature acting as a tool to examine the behavior of growth curves. Thus, 
any additive terms associated with geometrical structures is useful for adjusting a 
proper path describing the behavior of any epidemic curve. Since 1990s Antonelli's 
approach has been extending to include other types of geometries such as Finsler 
geometry rather than Riemannian [3]. Consequently, Antonelli and Bradbury (1997) 
have used Gompertz growth and allometric relationships to build a dynamical the- 
ory of ecology, evolution and development in colonial organisms using Finslerian 
Geometry. 

Accordingly, geometrization of Si-model can be considered as an introductory 
step to the geometric version of SIR-model and others. Moreover, this approach 
has been extended to include Finslerian Geometry due to its richness of imposing 
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many interacting factors geometrically rather than the Riemannian one [5]. The 
importance of geometrizing such a model is to construct the path equations and their 
corresponding path deviation equations to know the behavior of the growth curves 
and their effect after a slight perturbation on them. This can be seen by studying the 
stability conditions from examining the motion of path deviation equations. In the 
present work, we are going to utilize the nonlinear version of allometry, Riemannian 
and Finslerian, of Antonelli for obtaining the relevant equations of epidemics in 
Si-models using different types of geometries with some details on its extension to 
include SIR models and others by increasing the dimension of the manifold. This 
type of extending dimensions to examine various epidemic models will be examined 
in our future work. 



2 Historical Background 

2.1 SI-Model of Epidemics 

It is well known that the simplest model of describing population growth is stemmed 
from the famous prey- predator model cf[7] can be expressed as follows: 

and 

[1.6] 

whereS'(t) the susceptive class I(t), the infected class and a\ s a parameter. Equa- 
tions (l.a)and(l.b) can be related together using the following condition: 

S(t) + I(t) = N [2] 

where N stands for the total number of the population. The above model has been 
modified, Lotika-Volterra model , to accommodate some other factors based on some 
oscillations either internally or externally. In this study, it is sufficient to display 
the amended version owing to this extra factor to take the following form [8] 

< ^ = -a 1 SI + a 2 S [3.a] 

and ^ 

— = aiSI + a 2 I [3.6] 

where a 2 is another a parameter with constant coefficient, used to measure this 
interaction. From this perspective , it is vital to construct similar sets of path 
equations using Riemannian and Finslerian geometries respectively. 

2.2 The Concept of Allometry 

In 1936, Sir Julian Huxley introduced the concept of allometry, or the experimental 
study of relative growth of parts of animals, via log -log plots of measurements of 
morphological characteristics of individuals in a Euclidean geometry. These plots 
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resulted straight lines via statistical method of least squares. Around the same 
time Sir Joseph Needham made some experimental work to show that straight lines 
allometries. 

In 1944, J. Kittredge used the same concept to estimate the crown biomass of 
trees in forest stand by measuring the trunk girth, or diameter at breast height. 
Again some experiments of botanist J. Harper in 1962 who found Gompertz curves 
the best ones to fit for describing the growth of simple aquatic plants. 

By 1965 Laird studied on vertebrate growth using Huxlely's allometric law as 
well. Several applications of concept of allometry have been discussed in detail. 
Recently, Karl Niklas (1994) has recorded a variety growing plants and flowers 
satisfying the allometric concept as well [6]. 



3 Geometerization of Epidemics 

3.1 Underlying Geometry of Si-model: The Riemma- 
nian Version 

It is well known that, Antonelli and Voorhees (1974) have suggested the following 
metric in order to define geometrically the behavior of growth curves using allometric 
space [7]. 

g ij = e- 2a " xk 6 ij , [4] 
and its affine connection is given by 

r !jfc = -j9 l \9ik,\ + 9ik,l ~ 9jk,l), [5] 

where g % i is the the matrix inverse of gij. Suppose aj is a constant vector to define 
the Christoffel symbol in the following way: 

T l Ai = —on = constant [6. a] 

T]jj = aj = constant, i ^ j [6.6] 

I';,; "./,./ [6.c] 

F l j k = r 1 ^ = a.i = constanti ^ j ^ k [6.d] 
Thus, the metricity condition 

9ij;k = [7] 

becomes 

-2a m x m s:.. -2a m x m r y,nj -2a m x m c fii _ n 

e °Vk~ e ° mL k ~ e °nj L ki — U ' 

— za^Oij — e o m L k — e o n3 v ki — u, 

We can find out that the Riemann-Christoffel Curvature 

n abd — 1 ad,b - 1 ab,d + 1 ad L mb ~ 1 ab L md [ 6 - a J 
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has reduced to 



n abd — 1 ad L mb- 1 ab 1 md [°-°\ 

Thus the non vanishing components of the curvature tensor become 

R)il = OiiOiui^ j [9. a] 

MM 

R l jjt = ai aj,i ^ j ^ I. [8.c] 

3.2 Equations of Epidemics in Riemannian version of 
SI-Model 

It is well known that path equations of any growth curves can be obtained from the 
following the following applying variation principle on Lagrangian functions. This 
trend of imposing biological action principle has been stemmed prior to Euler by 
Mauphritis [2-7], 

L = g ab U a U b [10] 

However, in an approach to obtain path and path deviation equations from one 
single lagrangian using the Bazanski Lagrangian [10] which has the advantage that 
we obtain path and path deviation equations from the same Lagrangian: 

where a, b = 1, 2, 3, ..n. By taking the variation with respect to the deviation vector 
^ c and the tangent vector U c , and is covariant derivative with respect to a 
parameter s. From this perspective, it can easily be found that 

fJTTC 

— + T c ab U a U b = [12] 

and 

D 2 ^ c u A 

= R c abd U a U b * d [13] 

The path and path deviation equations of the Si-model can be obtained from 
the following Bazanski Lagrangian 

L = g»„U»— [14] 

By taking the variation with respect to the deviation vector ^ CT to get the fol- 
lowing path equation 

^ + Tl.S 2 + T l 22 I 2 + 2T\ 2 SI = 0, [15.a] 

and 

^ + T 2 n S 2 + T 2 2 I 2 + 2T 2 12 SI = 0, [15.6] 
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While taking the variation with respect to velocity vector U a to get their corre- 
sponding path deviation equations : 



and 



D 2 ^ 
~Dt 2 

D 2 m 



Dt 2 



s - R\ 12 SI^! + R\ 21 SW S + R\ 12 SI^i + BhiSIVs, [16.a] 
Rl 12 SI^! + R 2 121 SIV S + RlnSI^i + RhiSI^s, [16.6] 



i.e. 



and 



^f = (a 2 +2a 1 a 2 + a 2 )SI(^ s -^ I ) [16.c], 

= {a\ + 2 ai a 2 + a 2 2 )SI(^ s - */) [16.d] 

Consequently, we can suggest the path path and path deviation equations of a 
modified SI model can be following Lagrangian: 

L = g^^ + \W [17] 



to give, 



r1TT a 

dS 



+ T a b cU b U c = X a U a , [18.o] 



i.e 



and 



^ + r i 1 5 2 + 2ri 2 5/ + ri 2 / 2 = A5 [18.6] 



^- + T 2 11 S 2 + 2r 2 12 SI + T 2 22 I 2 = XI, [18.c] 
in which can be expressed as follows 

^ + 2/3SI - al 2 + aS 2 = XS [I8.d] 

and 

^ + 2(a - PL) SI - (P + aL)S + (a - /3L)/ 2 = A/. [18.e] 

And their corresponding path deviation equations become as follows: 

D 2 ^ a , A D^ a 

= R a bcd U b U^ d + A— [19.a] 

+ A ^ = ^12^/ + i2i2i5/*s + #212^*/ + RkiSWs, [19.6] 

and 

15^~ + A ^ = R ™ siiSl + R in si *s + #li 2 ^*/ + ^22i5/*s, [19.c] 
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which can be expressed as follows: 



+ X ^w = (a ' + 2aia2 + ^ SI ^ s - t 19 -^ 

and 

+ A ~d/ = («i + 2 «i«2 + o%)SI(* s ~ */) [19.e] 

3.3 Berwald type 

The metric tensor of Berwald type in Finslerian geometry can be described as follows 

9y = did j (^F 2 ) 

such a metric tensor is called a Minkowoski space with Finslerian norm ds = 
fidx 1 , dx 2 ). Thus, its corresponding geodesic can be described as 

^ + f} fc (x >y ) = [20] 

such that 

fjh = ^9 tm {Sh9jm + tfjflWi - S m 9jh), [21] 
where is a partial derivative with respect to non-linear connection. 

6 h = d h -r jh (x)^ K , [22] 

/* _ 1 n ii ( d 9ik ,dgu dg jk 

The Berwald type has some associated curvatures are defined in the following 
way [7]: 

B ljh = R ljh + G r j Dl rhk ~ G k D lhj P 4 ] 



where 



and 



such that 



to! dy ai dy a2 dy a3 <9y Q " 



x aia2a3....a m ™l a„.rvi £>„.rvo A„.rvs fl„.rv™ " ' 
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3.4 SI-Model in Finsler Geometry 

In 1991 Antonelli developed such a metric for 2-dimensional Berwald space with 
locally constant coefficients to become [3] 

F 2 = e 2a iX ^+l) (( ^1 )2 + { X 2 ) 2 ),i = l,2 [26] 

where X 1 and X 2 are Cartesian coordinates on R 2 and L acts as a perturbation 
parameter. The above relation can be related to Riemananian geometry by relaxing 
the term L 

- 9ij (X,X) = e 2 ^ L2 ^g i3 [27] 

The above system of equations together with their corresponding deviation vec- 
tor equations can be obtained from taking the action principle to the following 
Lagrangian. 

L BF = gij {x,y)X 1 -^- [28] 

Such that 

nA* fMr* ..... . „ . , 

_ _ _|_ f y\,J / •/• f •' iTf / 1 

where 



+ r; h w ft + c% [29] 



= ^. [30] 
ot 



Accordingly, geodesic equation may be as follows 

d 2 x i 



■ V' ll ,{.r. ! j) !J ' !l h .). 

and 



dt 2 

i.e. 



+ 2(/5 + qL),S/ + (/3L - q)/ 2 + (a - /3L),S 2 = [31.a] 



^ + 2(a - j3L)SI -{f3 + aL)S + (a - pL)I 2 = [31.6] 



and its deviation equation becomes 



B a bcd U b U c i>\ [32.a] 



Dt 2 

to become 



and 



i.e. 



Ds 2 



= B^SI*/ + ^2^7*5 + B^SI^i + B\ 21 Sli<s, [32.6] 
£ 2 12 S/*/ + B 2 21 tSW s + B 2 12 Sl4> T + Bf 21 SItf s , [32.c] 



1324 



7 



and 

-^r— - = 0. 32.e 
Ds 2 1 J 

Similarly, applying Antonelli's method to define Volterra's equation in ecology 

[2-7] , we obtain the its corresponding part of the Si-model in Berwald space, 

and 

d 2 I ,dS N o ~ ,dS.,dI. dl ,n ,dl , 

^ + r? l( -) 2 + 2r; 2 (-)(-) + rl 2 (-f = a(-), [33.6] 

Thus, the coefficients of its affine connection become: 

Y\ x = ai - a 2 L, [34.a] 

f^ 2 = -(ai-a 2 L), [34.6] 

r 12 = r 2i = -(«2 + ai^) [34.c] 

f? 2 = f 2i = -(ai -a 2 L) [34.d] 

ff 1 = -(a 2 + aiL) [34.e] 

fl 2 = -(a 2 + «!L). [34./] 
Thus, the components of path equations are expressed as follows: 

A Q 

— + 2{j3 + aL)SI + {pL - a){lf + (a - (3L){S) 2 = \S, [35.a] 

and 

^ + 2{a - QL)SI -{13 + aL)S 2 + {a - I3L)I 2 = XI. [35.6] 

Equations (35.a)and (35.b) are obtained by taking the variation with respect the 
corresponding deviation vector of the following modified Bazanski Lagrangian 

L = g a bX a —7— + \^ a X a [36]. 

Also, after some multiplications, we can find Their corresponding deviation equa- 
tions by taking the variation with respect to X c to become: 



Ds 2 Dt 

and 



+ A— f = Bl 12 SI*! + Bl 21 Sl4? s + Bl 12 Siif! + B^SIVs, [37.a] 



+ A— -i = 5 2 12 S/*/ + Bl 21 SIV s + BiuSI*! + B 2 2 21 5/^ s [37.6], 



i.e. 



Ds 2 £>t 

and 

L> 2 ^r 



' ' A^/ = [37.d], 



Ds 2 Dt 
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3.5 Antonelli-Finsler Metric Function 

In order to generalize the previous metric in (26) Antonelli has modified Finsler 
metric to become [7]: 

F = e*(J2 n i=1 (iT)™ , [38] 

where m > 2 and n = 2 . Thus their corresponding components of its affine 
connection are given as : 

rh=ai-^i(|)', [39.a] 

^i = ^i(|)3 + ^7i2, [39.6] 

rL = |/5i(f )«, [39.c] 

fl 2 = cr 2 -^ 2 (|)t. [39.d] 

Consequently, the Bazanski method to obtain path and path deviation equations 
will become in the following way: 
[i] The components of Path Equations: 



and 



at m — 1 to — 16 



§ + a 2 / 2 + + ^(ir- 2 / 2 = A 2 I. [40.a] 

at m — 1 to — 15 



[ii] The components of Path deviation Equations: 



+ A— - = 0, [41.o] 
Ds 2 Dt ' 1 J 

and 

It is well known that the interaction between 5&7 can be vitally important when 
to > 2 and the cases of increasing dimensions such as the following section to 
examine the possibility to geometrize SIR models. 



4 Geomertization of SIR Model 

In a similar way, we can extend our study to examine SIR model using the Bazanski 
method in each Ricmannian and Finslerian to become: 
(i)For Riemannian Geometry 

^ = [42.„] 
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and 

^ = 0, [42,, 

provided that 

S + R + / = N [43] 

where R is the number of recovered or dead group. 

The above equation can easily be obtained by assuming the following Lagrangian: 

L = g^^ + a {p) ^U v , a, b = 1,2, 3 [44] 
where is an arbitrary constant. 

From this perspective, we can develop the SIR model in its geometric version as 
follows 

~dT = a M u " ^ 

and their corresponding deviation equations become: 

-^ +a ^^r = R ^ UUUPq,a - ^ 

(ii) For Finslerian Geometry (Berwald Type): 

In a similar way to equations [42a-45] , we can easily construct the following 
equations: 

2^ = 0, [47.a] 
Dt 

^- = al, [47.6] 
Dt 

and 

^ = 0. [47.c] 
Dt 

Equations [47 a-c] are obtained by taking the action of the following Largangian: 

- bi> v 

L = g^^ + a M W, n,v= 1,2, 3. [48] 
We also can find their corresponding deviation equations to become: 

Dt 2 Dt 1 J 

For a complete description of this model will be examined in our future work. 
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5 Discussion and Concluding Remarks 



The paper deals with geometrization some epidemic models using their correspond- 
ing path equations. Also, we have obtained their corresponding path deviation 
equations from one single Lagrangian for different types of geometries based under 
a positive definite metric with constant affine connections and negative curvature. 
This process is based on a geometrized vision of concept of allometry to examine 
epidemic curves as by finding their path equations. The paper has also opened the 
window to impose non conventional types of geometries in future work to describe 
more complex models of SI or briefly in SIR model. One of good results is increasing 
the dimensions with preserving the concept of a positive metric with constant affine 
connections and negative curvature will include many other epidemic models such 
as susceptive-exposed- infective-recovered models (SEIR). 

Finally, this study can be extended to obtain such an appropriate prospective 
vision to visualize how an epidemic behaves on a longer times without many pa- 
rameters as needed in case of the traditional regression analysis. This is not the 
optimal case to get an exact form of epidemic curve, but it is a tend to apply such 
a concept of geometrization. In addition in our future work we may in need to de- 
velop these geometries to include some new parameters affecting the Si-model such 
as transmission rate [12] to become: 

^ = -aS p P [50] 
dt 1 J 

can be studied in future to be geometrized. 

Also, for the SIR model can also be extended to include some interactions like 
population fertility or immigration to be considered in demographic using partial 
differential equations of McKendrick- von Forester with boundary conditions [13]. 

d d 

(t^ + T^)n{a,t) = -fin(a,t) + I{a,t), 0<a<u,t>0 [5] 



such that 

n 



(0,t) = B(t)= m(a)n(a,t)da, t>0 
Jo 

n(a, t) = no(a), < a < to 



where uj represents the maximum life span of individuals, B(t) is the number of 
newborn per unit time t , m(a) is a mortality rates and fi(a) is a constant. The 
evolution of density n(a, t) of individuals aged a and timet, will open the window 
to be included geometrically as non linear connections in Finslerian approach. 
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